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5_i ■ Abstract 



There is a well known anomaly between the value of the Fermi decay constant 
extracted from super-allowed Fermi beta-decay of nuclear isotriplets and that re- 
quired by unitarity of the Cabibbo-Kobayashi-Maskawa matrix. This discrepancy 
remains at the level of a few tenths of a percent after the most rigorous investiga- 
tion of conventional nuclear and radiative corrections. Within the framework of the 
quark-meson coupling model of nuclear matter, which has been previously applied 
successfully to phenomena such as nuclear saturation and nuclear charge symmetry 
violation, we show that it is possible to understand a significant fraction of the 
observed anomaly. 
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It is clearly very important to refine our understanding of the weak coupling to quarks 
as much as possible. Testing the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix is one of the more challenging aspects of this general problem. In particular, 
the most accurate experimental measurement of the vector coupling constant in nuclear 
beta-decay comes from super-allowed + -0 + transitions between nuclear isotriplet states. 
However, in order to relate these precise measurements to the quark-level vector coupling, 
V u d, one needs to apply a number of small nuclear structure corrections in addition to 
the relatively standard radiative corrections 0. Despite intensive study of these nuclear 
"mismatch" corrections |3], |4], ^ |6| there remains a systematic difference of a few tenths 
of a percent between the value of V u d inferred from the vector coupling measured in muon 
decay, G^, and unitarity of the CKM matrix and those determined from the nuclear ft- 
values. For recent summaries we refer to the reviews of Wilkinson and Towner and 
Hardy ||, and also to the recent report by Savard et al. || of accurate data on 10 C. 

Until now the nuclear corrections have been explored within the framework of con- 
ventional nuclear theory with point-like nucleons. Of course, for the nucleon itself there 
has been considerable investigation of the effect on the vector form-factor of the breaking 
of CVC caused by the small u-d mass difference in QCD. The Ademollo-Gatto theo- 
rem 



ID, |TT| tells us that any corrections must be of second order in (m^ — m u ) - a result 



that has survived [|P2| , |T3|1 suggestions that it might fail because of p-u mixing [|T4| . While 
this is necessarily very small, the measurements of V u d and are also extremely precise. 
Thus we have been led to ask whether this small nuclear discrepancy might be associated 
with a change in the degree of non-conservation of the vector current caused by nuclear 
binding. 

In order to investigate whether nuclear binding might influence the Fermi decay con- 
stant of the nucleon itself one needs a model of nuclear structure involving explicit quark 
degrees of freedom which nevertheless provides an acceptable description of nuclear bind- 
ing and saturation. The quark-meson coupling (QMC) model of Guichon [15|] seems 
ideally suited to the problem. In this model, nuclear matter consists of non-overlapping 
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nucleoli bags bound by the self-consistent exchange of a and uj mesons in the mean-field 
approximation. It has been extended by Yazaki et al. to include a centre of mass 
correction and by the present authors to include the p and an isovector-scalar meson (the 
5) [|T7| , 18] . As well as providing an excellent description of the properties of nuclear mat- 



ter, it has been applied successfully to the calculation of nuclear structure functions [IS 



and (allowing for a difference in the quark masses, m„ 7^ m^) to the problem of the 
Okamoto-Nolen-Schiffer anomaly |2(| in mirror nuclei jrj|. Furthermore, the relationship 
between the QMC model and Quantum Hadrodynamics has been investigated. The 
fascinating result is that for infinite nuclear matter the two approaches can be written 
in an identical form, except for the appearance of the quark-scalar density in the self- 
consistency condition for the scalar field |17], p2[] . The simplicity of this finding suggests 
that it may be rather more general than the specific model within which it was derived. 

We shall briefly review the main features of the model with unequal quark masses 
before turning to the specific problem of the Fermi form-factor. Suppose that the mean- 
field values for the a, u (the time component), p (the time component in the third direction 
of isospin) and 5 (in the third direction of isospin) fields, in uniformly distributed nuclear 
matter with N Z, are a, uj, b and 5, respectively. The nucleon is described by the static 
spherical MIT bag in which quarks interact (self-consistently) with those mean fields. 
Then the Dirac equation for a quark field, ij), in a bag is given by 

[17 • 9 - (mi -V a - ^T Z V 5 ) - 7°(K, + ^TzVpM/j = 0, (1) 

where V a = g%cf, V w = g^u, V p = g^b and V$ = gjS with the quark-meson coupling 
constants, g$, g%, g q p and g q s . The subscripts, i and j, denote the i-th quark (i=u or d) in 
the proton or neutron (j=p or n). 

The normalized, ground state wave function for a quark in the nucleon is then given 

by 



3 



where 



6i/j = Qi/j + Rj(V w ± ~V P ), for a rj quark (3) 
tfjf = 2R 3 j f (x i )[Q i/j (Q i/j - 1) + i^m*/2]/a?, (4) 
A/j = ^ (fij/j - Rjm*)/ (tti/j + Rjmt), (5) 



with Vti/j = yx? + (Rjin*) 2 and Xi the quark spinor. The effective quark mass, m*, is 
defined by 

m* = nii- (V a ± for a ^ quark. (6) 

The linear boundary condition at the bag surface determines the eigenvalue Xi. 

Taking the spin-flavor wave function for the nucleon to have the usual SU(6) form, 
the nucleon energy is given by E\ ag + 3V U ± ~V P for PJ, where the bag energy is 



EL, = E '"^" 2 ° + |xBi§, (7) 



B is the bag constant and Zq is a phenomenological parameter initially introduced to 
account for zero-point motion. Here we use it to correct for spurious cm. motion as well 



- rather than following the more elaborate procedure of Ref. JT6[. The reason is that, as 
illustrated by the Ademollo-Gatto theorem, the deviation of the vector form-factor from 
unity is a very subtle effect. This method is the only one in which we have been able to 
guarantee that in the free nucleon case the deviation is proportional to (rrid — rn u ) 2 . The 
effective nucleon mass, M*, in nuclear matter is given by minimizing eq.(^) with respect 
to Rj. 

To see the sensitivity of our results to the bag radius of the free nucleon, we choose 
m u = 5 MeV and vary the parameters, B, zq and m^, to fit the physical proton and 
neutron masses for several values of the average, free bag radius i? (= 0.6, 0.8, 1.0 fm). 
Since the electromagnetic (EM) self-energies for p and n contribute to the masses we 
adjust the parameters to fit the bare proton mass, M p = 938.272 — 0.63 MeV, and the 
bare neutron mass, M n = 939.566 + 0.13 MeV, where +0.63 MeV and —0.13 MeV are the 
EM self-energies for p and n, respectively ||23|| . We then find that B 1 ' 4 = (210.7, 169.8, 
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143.6) MeV, z = (4.011, 3.305, 2.600), and m d = (9.252, 9.242, 9.233) MeV for R = 
(0.6, 0.8, 1.0) fm, respectively. It is crucial for the present application that the bag radius 
of the proton is a little smaller than that of the neutron in matter and that this difference 
increases slightly with density. 

At the present time the extension of the QMC model to finite nuclei is not com- 
plete [ 24f . In particular, the variation of the scalar and vector fields across the finite size 



of the nucleon presents quite a challenge. Thus, for the present we are constrained to 
discuss the problem within the framework of infinite nuclear matter. In this case we take 
the Fermi momenta for protons and neutrons to be kp p and kp n , respectively. These are 
defined by p p = kp p /(3n 2 ) and p n = kp n /(3n 2 ), where p p and p n are the densities of p 
and n, respectively, and the total baryon density, pb, is given by p p + p n . The u field 
is determined by baryon number conservation, and the p mean-field by the difference in 
proton and neutron densities (ps below). On the other hand, the scalar mean-fields, a 
and 5, are given by self-consistency conditions (SCCs) [[H], Since the p field value is 



given by b = g p p 3 /(2m 2 ), where g p = g q and p 3 = p p — p n , the total energy per nucleon, 
E to t, can be written 



Etot = E P dkJMf + P + ^a 2 + ^±6 2 + ^p B + Jf-,,,. , * , 

Pb(2tt) 3 .±£ n J V 2p B 2p B 2ml 8m 2 p p B 

where g w = 3g«. 

We determine the coupling constants, g 2 /^ = (22.62, 20.78, 19.47), where g a = 3g%, 
and gl/An = (5.547, 4.637, 3.979) for R = (0.6, 0.8, 1.0) fm, respectively, so as to fit the 
binding energy (—16 MeV) and the saturation density (po = 0.17 fm -3 ) for equilibrium 
nuclear matter. It is important to stress that the S plays an insignificant role in our 
present calculation, while the p has no influence in symmetric matter. Nevertheless, for 
the record we note that g 2 /4:ir is chosen to be 2.82 |25| (where gg = gf) and the p meson 
coupling constant, ^/4vr, is chosen to be (5.223, 5.340, 5.423) for R = (0.6, 0.8, 1.0) fm, 
respectively, so as to reproduce the bulk symmetry energy of nuclear matter, 33.2 MeV. 
One of the important consequences of including the internal structure of the nucleon is 
that the model gives a very respectable value for the nuclear incompressibility, around 
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Figure 1: I U ci(pb) / 'Aid(O) as a function of the nuclear density. The solid, dotted and dashed 
curves are for Rq = 0.6, 0.8 and 1.0 fm, respectively. 

200 ~ 300 MeV. 

Armed with an acceptable, quark-based model of nuclear matter we can now inves- 
tigate the variation with density of the quark vector current matrix element. For Fermi 
beta-decay we need the matrix elements 

Iw(Pb) = / dVipl /p ipi,/ n , (9) 

J Bag 

with i' = d and i = u for the d — >• u conversion and i = i' = u or d for the two spectator 
quarks. As the radius of the proton and neutron are different we integrate over the 
common volume. 

As remarked earlier, the first test of the consistency of the calculation is that at zero 
density the deviation of the vector form-factor from unity is quadratic in — m u . By a 
numerical check we can see that this is indeed the case . 

In fig.|l| we show the ratio I u d(pB)l I u d{fy, for our three choices of the free (average) bag 
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radius (in symmetric nuclear matter), as a function of ps- We also give full details of those 
values at several nuclear densities in table |l|. The decrease in In* as the density increases 
is a direct consequence of the increasing difference between the proton and neutron radii 
- that is the smaller volume of overlap. From the figure and the table, as the density 



Table 1: Variations with density of Iw for three bag radii. 

.Ro(fm) Pb/po lud hu hd 

0.999993 0.999995 0.999995 

0.6 0.5 0.999869 0.999870 0.999871 

1.0 0.999751 0.999751 0.999752 

0.999991 0.999994 0.999994 

0.8 0.5 0.999840 0.999842 0.999843 

1.0 0.999702 0.999703 0.999705 

0.999987 0.999992 0.999992 

1.0 0.5 0.999814 0.999817 0.999819 

1.0 0.999664 0.999666 0.999669 



increases from zero, we see that the deviation of Iw(pb)/ Iu'{0) from unity is roughly 
linear in (The deviation is also roughly linear in irtd — m u .) We can then summarise 
the results as: 

^Ui- «.x(<*y do) 

hi'{0) \poJ 
with aw ~ (2.4,2.9,3.3) x 10~ 4 for R = (0.6, 0.8, 1.0) fm, respectively (for any combi- 
nation of ii') . 

The evaluation of /t-values involves the inverse of the product of I u d, I U u and Idd 
squared. Since for a given, free (average) radius of the bag each of these matrix elements 
decreases by roughly the same amount, the fractional increase in the ft- value with density 
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is therefore 

with 6 approximately six times the decrease in each integral - i.e. b ~ (1.5, 1.8, 2.0) x 10~ 3 
for Rq = (0.6, 0.8, 1.0) fm. Thus the increase in the /i-value at po/2 ranges from 0.075% 
to 0.10%, while at po it lies between 0.15% and 0.20%. This is to be compared with a 
violation of unitarity of the CKM matrix of 0.35 ± 0.15% in the most recent analysis of 
Towner and Hardy H. 

It is not possible to draw unambiguous conclusions from a comparison of theoretical 
results in infinite nuclear matter with data from finite nuclei. At po/2 the calculation 
suggests a reduction in the violation of unitarity by about 1/3, while at po a correction as 
big as 0.2% brings the discrepancy back to only one standard deviation. We also note that 
this result is rather sensitive to the proton fraction, f p = p p /pb- The overlap integrals 
with f p less (greater) than 0.5 decrease more slowly (rapidly) than in the symmetric case. 

Before concluding our discussion it might be helpful to make a qualitative observation 
concerning the validity of our results. The essential physics involved in our calculation 
is charge symmetry violation, in particular, the fact that in nuclear matter the confining 
potential felt by a quark in a proton is not the same as that felt by a quark in a neutron. 
We have already explained that a relativistic field theory only yields the right order of 
magnitude for nuclear charge symmetry breaking if the relevant mass scale involves quarks 
rather than nucleons f26|. In this sense the ONS anomaly may prove to be something of 
a "smoking gun" for quark degrees of freedom in nuclei. This is even more obvious here; 
it is only because the nuclear charge symmetry violation occurs at the quark level that 
it can produce a deviation of the vector form factor of the bound nucleon from its free 
value. 

In summary, within a specific, relativistic quark-meson coupling model of nuclear mat- 
ter we have found that the change of the internal structure of the nucleon itself produces 
a correction to the Fermi beta-decay form-factor of the sign and order of magnitude re- 
quired for consistency with the unitarity of the Cabibbo-Kobayashi-Maskawa matrix. It 
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is important that the model dependence of this result be further examined as soon as 
possible - although we have given arguments why the qualitative result may be quite 
general. Finally we note that the present model deals only with nuclear matter and it is 
vital that the formalism be extended so that we can apply the model quantitatively in 
finite nuclei. 

We would like to thank Peter Jackson for raising this issue with one of us during the 
celebration of Erich Vogt's 65th birthday at TRIUMF last December. We would also like 
to acknowledge helpful discussions with T. Hatsuda, I. Towner and D. Wilkinson during 
WEIN'95. This work was supported by the Australian Research Council. 
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